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ABSTRACT 

When a magnetically dominated super-fast-magnctosonic long/soft gamma-ray burst 
(GRB) jet leaves the progenitor star, the external pressure support will drop and the 
jet may enter the regime of ballistic expansion, during which additional magnetic accel- 
eration becomes ineffective. However, recent numerical simulations by Tchekhovskoy 
et al. have suggested that the transition to this regime is accompanied by a spurt of 
acceleration. We confirm this finding numerically and attribute the acceleration to a 
sideways expansion of the jet, associated with a strong magnctosonic rarefaction wave 
that is driven into the jet when it loses pressure support, which induces a conversion 
of magnetic energy into kinetic energy of bulk motion. This mechanism, which we dub 
rarefaction acceleration, can only operate in a relativistic outflow because in this case 
the total energy can still be dominated by the magnetic component even in the super- 
fast-magnctosonic regime. We analyse this process using the equations of relativistic 
MHD and demonstrate that it is more efficient at converting internal energy into ki- 
netic energy when the flow is magnetized than in a purely hydrodynamic outflow, as 
was found numerically by Mizuno et al. We show that, just as in the case of the mag- 
netic acceleration of a collimating jet that is confined by an external pressure distribu- 
tion - the collimation acceleration mechanism - the rarefaction-acceleration process in 
a magnetized jet is a consequence of the fact that the separation between neighbouring 
magnetic flux surfaces increases faster than their cylindrical radius. However, whereas 
in the case of effective collimation-acceleration the product of the jet opening angle 
and its Lorentz factor does not exceed ~ 1, the addition of the rarefaction- acceleration 
mechanism makes it possible for this product to become 3> 1, in agreement with the 
inference from late-time panchromatic breaks in the afterglow light curves of long/soft 
GRBs. 
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1 INTRODUCTION 

In the 'standard' model of lo ns-duration, soft-spectrum 
gamma-ray bursts (GRBs; e.g. IPiranl l2005h . the prompt 
high-energy emission arises in ultra-relativistic (bulk 
Lorentz factor F>10 2 ), highly collimated (opening half- 
angle of a few degrees) jets. The high Lorentz factors 
are inferred from the requirement of a sufficiently low 
opacity to photon-photon annihilation or to scattering by 
photon annihilation-p roduced electron-positron pairs (e.g. 
iLithwick fc Sarill200ll ). The high collimation makes it possi- 
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ble to reduce the total flow energy down to values that are 
comparable to the energy of stellar explosions. 

As the jets make their way through the interstellar 
medium they sweep the ambient gas into a thin relativistic 
shell that itself becomes a strong source of electromagnetic 
radiation. After a sufficiently large amount of gas is swept 
up, the shell begins to decelerate. Simplified models of this 
interaction predict that a panchromatic break would then 
occur in the afterglow light curve, provided that TOj 3> f , 
where Oj is the jet opening half-angle. The detection of such 
'jet breaks' on time scales of > 1 day in several (primar- 
ily long/soft) GRB afterglows has been a key reason for 
the widespread adoption of the jet model for these sources. 
The observed break parameters have made it possible to 
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Figure 1. Lorentz factor and magnetic field lines in models Bl (left panel) and Bib (right panel). 



deduce the opening half-angle of these jets. The inferred 
values are strongly model depend ent, although they usu- 
ally l ie the range of 1° - 10° (e.g. lRhoads|[l999l : ISari et al.1 
Il999l ). Although more recent observations by the Swift satel- 
lite have revealed that late-time panchromatic jet breaks are 
not that common and that vari ous aspects of the jet model 
may need to be modifi e d (e.g. Meszarosl 20061 ; IPanaitescul 
120071 : iLiang et al.1 |200S| ; iRacusin et al.l l2009h . these breaks 
remain the strongest evidence for collimated outflows in 
GRB sources. 

The supernova connection of long/soft GRBs provides 
strong support for theoretical models of their jet engines 
that invoke dying massive and rapidly rotating stars. These 
models are generally divided into two groups depending 
on the mechanism of jet acceleration. In the first group 
the acceleration is driven by the thermodynamic pressure 
of plasma heated to ultra-relativistic temperatures via the 
annihilation of neutrinos and antineutrinos emitted by the 
accretion disk formed around the central black hole, thus 
tapping the disk thermal energy. In the second group the 
acceleration is driven by magnetic stresses, tapping the ro- 
tational energy of the disk or the central compact object 
(neutron star or black hole). At present, both the magnetic 
and the thermal mechanisms seem possible, although the 
lack of detection of a thermal component in the spectra of 
some GRBs is consistent with the notion that, at least in 
certain cases, the outfl ow is initially magnetically dominated 
ijZhang fc Pe'erll2009l ). 

The magnetic jet acceleration mechanism has been 
the subject of theoretical study for many years. Due 
to the mathematical complexity of magnetohydrody- 
namics (MHD), which is even more pronounced in the 
relativistic limit, it has been possible to find analytical 
and semi-analytical solutions only for a rather limited 
number of problems characterized by a high degree of 



symmetry. In fact, there is only one available exoc£] 
solution of the relativistic MHD equations including 
thermal and magn e tic eff ects, the self-similar model of 
IVlahakis fc Konigj (|2003l ). The advance of numerical 
methods for relativistic MHD during the last decade has 



significant progress (e.g. Komissarov 2001 


. 2004; McKinnev 


200£ 


; komissarov & Barkov 120071; 


Komissarov et al. 


2007 


; Tchekhovskov et al.l |2008|; iBucciantini et. al. 


2008 


; Komissarov et al.l 20091; Tchekhovskov et al 


2009a; 


Komissarov & Barkovl |2009|; Tchekhovskov et al.l 


2009b; 


McKinnev & Blandford 20091; Bucciantini et. al.ll2009T). 


In particular. Komissarov et al.l ( 2009h investig 


ated the 



magnetic acceleration of ultra-relativistic flows within chan- 
nels of prescribed geometry, 2 oc r a (where r and z are 
dimensionless cylindrical coordinates) , determined by the 
shape of coordinate surfaces of elliptical coordinates. Such 
shapes correspond to power-law distributions of the confin- 
ing pressure that can approximate the expected distribu- 
tions in the envelopes of the progenitor stars in both the 
collapsar and the magnetar scenarios (which correspond to 
the dying star leaving behind a black hole or a neutron star, 
respectively) as well as the effect of a confining disk wind. 
Among other results, they found that in the case of a gradu- 
ally widening channel, a < 1, the acceleration is not efficient, 
whereas in the case of a channel with gradually increas- 
ing collimation, a > 1, the acceleration is effective but the 
asymptotic flow obeys the inequality F9j < 1. These numer- 
ical results have been strengthened by complementary the - 
oretical analyses (|Komissarov et al.ll2009l ; lLvubarskvlf2009l ) , 
which led to the question of whether the magnetic models 



1 Exact in the sense that the full system of MHD equations - 
including the fluid-inertia terms and the trans-field component of 
the momentum equation - is integrated. 
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Figure 2. Left panel: Evolution of V, oT = fj, — V and fi along the magnetic surface enclosing 80% of the total magnetic flux = 15) 
in models Bl (solid line), Bla (dash-dotted line), and Bib (dashed line), with the distance measured in units of the light-cylinder radius 
n c of the source. Right panel: Evolution of T9 V , where 8 V is the opening half-angle, along the jet in model Bib. The solid lines show the 
variation of this parameter across the jet as a function of the poloidal magnetic flux >!/ for z = 1, 2, 4 and 8 X 10 r\ ci with higher values 
of z corresponding to higher curves at \P = 15. The dashed curve shows V8 V at z = 8 X 10 5 r\ c in model Bl. 



can accommodate the jet breaks at a ll. Recently, a simi- 
lar num erical study was carried out by Tchekhovskov ct al 
j), who confirmed the findings of iKomissarov et al 



They have, however, also considered a somewhat 
modified setup, wherein at some distance from the origin, 
roughly corresponding to the stellar surface, the channel 
geometry changes from progressively collimating to pro- 
;ressive ly de-collimating. In this setup rTchekhovskov et akl 
2009bl ) observed a remarkable change in the jet behaviour: 
beyond the point where the geometry changed, the jet speed 
underwent a strong boost that was accompanied only by a 
very small increase in the jet opening angle. As a result, 
the asymptotic flow had T8j 3> 1, which made it possible for 
late-time jet breaks to occur. The exact shape of the channel 
above the transition point did not seem to matter. 

In this paper we describe (Section [2| simula tions that 
confirm the results of lTchekhovskov et al. (2009b]) and then 
analyse (Section [3]) the underlying physical mechanism. We 
discuss and summarize our results in Section [4] 



2 NUMERICAL SIMULATIONS 

The numerical metho d is exactly the same as in 
IKomissarov et al. (2009) and we refer the reader interested 
in technical details to that paper. Here we only remark that 
our numerical code is based on the Go dunov-type scheme for 
relativistic MHD (|Komissarovl 1 199*91 ) and that we look for 
steady-state axisymmetric solutions using time-dependent 
simulations with time-independent boundary conditions. 

For our purpose we selected model Bl from 
IKomissarov et all (|2009i ), for which the channel shape pa- 
rameter is a — 3/2. This model describes a cold flow with 
values of the field-line constant /i, defined as the energy flux 
per unit rest-mass energy flux, as large as fj, = 620. This 
constant sets the upper limit on the Lorentz factor that can 



be achieved in this model via ideal MHD mechanism0 The 
initial ratio of the Poynting flux to the hydrodynamic en- 
ergy flux is (Jo = jU/To — 1 ~ /i, as the initial Lorentz factor, 
To, is close to unity, corresponding to a sub-Alfvenic flow. 
The base rotation is uniform, with the dimensionless light- 
cylinder radius n c being ~ 1.6; we used the distance of the 
inlet boundary from the origin as the unit of length. 

To examine the effect of changing the channel shape, 
we map the solution at z ~ 10 3 and 7 x 10 4 onto the inlet 
boundary of a new grid corresponding to a conical channel 
of the same local radius and with a vertex located at the 
origin. In this way we introduce a change in the channel 
shape to a profile that results in less collimation. We then 
proceed with the simulations on the new grid, following the 
same procedure as in model Bl. The two solutions obtained 
in this way, which we denote as Bla and Bib, are analysed 
below. 

Fig.rTJshows the overall geometry as well as the shape of 
the magnetic surfaces and the evolution of the Lorentz factor 
in models Bl and Bib. One can see that, in contrast with the 
situation in model Bl, the field lines in model Bib straighten 
out. One might think that this reflects the conical shape of 
the channel, but this is not so. In fact, the jet is separated 
from the channel wall by a near vacuum0 This is the reason 
why the red-coloured boundary layer in the right panel of 
Fig. [1] is free from magnetic field lines. A simila r sepa ration 
has been seen in model E of IKomissarov et al.l (|2009l ) and, 



2 In terms of the energy budget, magnetic acceleration represents 
the conversion of Poynting flux into kinetic energy. The ratio of 
the kinetic energy flux to the rest-mass energy flux is equal to 
the bulk Lorentz factor of the flow. Thus, when the Poynting flux 
in the jet is fully converted into kinetic energy, the acceleration 
process stops and one has fi = F. 

3 In the simulations, the mass density in the vacuum zone is kept 
above zero by numerical effects. 
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we beli eve, also in the simulations of iTchekhovskov et al.l 
( |2009bh . As a result, the wall and the jet are causally dis- 
connected and the precise shape of the wall does not matter. 

In model Bib the jet Lorentz factor approaches its max- 
imum possible value T max = [i at the jet boundary, signalling 
a total conversion of the Poynting flux. The acceleration 
is weaker in the jet interior but, as one can clearly see in 
Fig- El it is still more effective than in model Bl. This is fur- 
ther illustrated in the left panel of Fig. [5] which compares 
the acceleration in both models along the magnetic surface 
enclosing 80% of the jet's total poloidal magnetic flux. In 
the case of model Bla, where the channel opens up much 
earlier, the jet also passes through a phase of rapid accel- 
eration. However, the acceleration slows down dramatically 
when the jet enters the phase of ballistic expansion. As a 
result, the final Lorentz factor in this model is even lower 
than in model Bl. As one can see from Fig. [T] the opening 
angle of the jet does not change much after passing the point 
where the channel widens, i n agreement with the results of 
ITchekhovskov et all (|2009bl ). Thus, the product F0 V , where 
6 V is the local opening half-angle, is expected to increase 
following the rapid increase of the Lorentz factor. This is 
indeed the case, as one can see in the right panel of Fig. [2] 
In model Bib this product is much larger than in model Bl, 
approaching values that are 3> 1 near the jet boundary. Con- 
cluding this section, we reiterate that o ur results are in very 
pod ag reement with those obtained in ITchekhovskov et al.l 
2009bh , thus confirming that the effect is real and must 
have a robust physical basis. 



3 ACCELERATION MECHANISM 

Once the jet separates from the wall in models Bla and 
Bib it enters the phase of free expansion and eventually 
becomes a ballistic conical outflow with radial streamlines. 
During the transition from the one regime to the other a 
strong fast-magnetosonic rarefaction wave propagates into 
the jet. Since at this point the jet is causally connected 

the wave is not confined to 



let, since at tms point 
(see iKomissarov et all 120091 ') , 
a boundary layer but propagates all the way to the jet core. 
This is clearly seen in Fig. which compares the distribu- 
tions of the magnetic pressure b 2 /8tt (where b is the magnetic 
field amplitude in the comoving frame) in models Bl and 
Bib downstream from the channel- widening point. The ob- 
served jet acceleration is apparently related to the properties 
of this wave. In fact, this phenomenon has already been seen 
in other numerical sim ulations of both purely hydrodynamic 
and mag netized flows (|Alov fc Rezzollall2006l ; | Mizuno et al.l 
l200ct ). Here we elucidate its physical nature by considering 
a simpler one-dimensional problem of relativistic expansion 
into vacuum in a slab geometry. The rarefaction wave in this 
case is described by a self-similar solution known as a 'sim- 
ple wave'. Although such a flow is not identical to the one in 
our simulated jets, it nevertheless captures the underlying 
physical mechanism. 



3.1 Simple waves 

Consider a quasilinear hyperbolic system 

(1) 



°2- 
o 




CO 



Figure 3. Propagation of a rarefaction wave across the jet in 
model Bib. The solid lines show the magnetic pressure dis- 
tribution as a function of the poloidal magnetic flux at z = 
1, 3, 9, 30 X 10 7"i c for model Bl, whereas the dashed lines show 
the corresponding distributions for model Bib . In this plot, the 
higher the value of z, the lower the corresponding curve. 



where P — (Pi, P2, . . . , P n ) T is a vector of dependent vari- 
ables, e.g. the gas pressure, density, etc., and A(P) is an 
(n x 71) matrix. When the initial (t — 0) configuration de- 
scribes two uniform states separated by a discontiniuty at 
x — (a Riemann problem), the system allows self-similar 
solutions that depend only on £ = x/t. In general, the solu- 
tion involves a combination of shock and rarefaction waves, 
but in the case of expansion into vacuum it describes a sin- 
gle rarefaction wave (a simple wave). In terms of the new 
independent variable £, the system of equations ([TJ reduces 
to 



(A- 



(2) 



where I is the unit matrix. This equation has the form of 
the characteristic equation for the matrix A, and thus £ is 
equal to one of the eigenvalues of A, 



C = HP) , 



(3) 



r(P) 



(4) 



with 
dP _ 

*T = 

being the corresponding right eigenvector. Each eigenvalue 
represent the phase speed of a particular wave, whereas 
the right eigenvector determines the connection between the 
variations of the dependent variables that are induced by the 
wave. Equation @ can be conveniently written as 



dPi_ 

ri 



dP 2 

T2 



dP n 



{->) 



dP dP 

dt dx 



Integrating this system, one finds Pi = Pi(Pi), where i — 
2 . . . n, and then equation ((3]> (the equation of characteris- 
tics) can be used to obtain Pi = 



© 0000 RAS, MNRAS 000, 000-000 



Rarefaction acceleration of magnetized GRB jets 5 



The systems of equations for relativistic hydrodynam- 
ics and for relativistic MHD can be written in the form of 
equation JTJ ■ The eigenvalues and eigenvecto rs for all hy - 
per bolic waves can be found, for example, in lAnilei (|l989T ) 
and Komissarovl ( 199£ ) . Simple waves have been analysed in 
iMartf fc Miiller ( 1994 ) for purely hydrodynamic flows and 



Romero et. al 



( 2005h for magnetized flows in which the 



magnetic field is tangential to the discontinuity and orthog- 
onal to the flow velocity. 

3.1.1 Case 1 

Consider a simple one-dimensional problem with plane ge- 
ometry describing the evolution of an initial discontinuity 
that separates a uniform magnetized cold plasma at rest on 
the left and vacuum on the right. This problem is related 
to the transverse expansion of our simulated jets after they 
lose external support, as seen in the jet comoving frame. The 
discontinuity decays into a single wave, namely the simple 
fast-magnetosonic rarefaction wave. The left front of this 
wave (the head) propagates into the left state with the lo- 
cal fast-magnetosonic speed and the right front (the tail) 
moves into vacuum with some finite terminal speed. In such 
a simple geometry, the Lorentz force only involves the mag- 
netic pressure, and the equations of MHD reduce to those 
of hydrodynamics of an ideal gas with a polytropic index 
7 = 2. In this case the effective sound speed a is equal to 
the fast-magnetosonic speed, which satisfies 



(!)' 



Arvpc 2 + b 2 1 + a ' ^ 

where p is the rest-mass density and a = b 2 /4-Kpc 2 is the 
local magnetization parameter. This property allows us to 
utilize the results obtained for the corresponding problem 
in relativistic gas dynamics. In particular, the integration of 
equation ([5| yields the integral across the rarefaction 



J. 



1 + v x /c\ fl + a/ 



(7) 

Vx/c J \ 1 — a/c J 

ijMartf fc Mtillerll 19941 ). where v x is the local flow speed. This 
equation can be used to deduce vo, the wave's expansion 
speed into vacuum. Specifically, by equating the value of 
the above invariant as derived from the given left state to 
the value obtained at the boundary with the vacuum (where 
a = 0), one finds 



«o = 



J+ + 1 



(8) 



r 2/ 2n-i/2 l + a 2 /c 2 

r = (l-«o/c) =- j-r^. (9) 

1 — a z jcr 

For example, if in the left state, where v x = 0, we have 
b 2 = 4-Kpc 2 , and hence a = c/-/2, then the Lorentz factor 
of expansion into vacuum is To = 3. Next, consider another 
reference frame that moves along the initial discontinuity 
with a Lorentz factor Tj - this corresponds to the jet source 
frame in our simulations. The Lorentz factor of plasma at 
the right front of the rarefaction wave in this frame is 



r = r,-r 



(10) 



which in our example is three times larger than Tj. Thus, 
the seemingly weak acceleration in the jet frame may corre- 



spond to a huge boost in the jet source frame: in our example 
the Lorentz factor increases from Tj = 200 to the left of the 
rarefaction wave to T = 600 at the boundary with vacuum. 
This is the essence of the acceleration mechanism that op- 
erates in the free expansion regime of the jet simulations. 
However, the jet expansion problem cannot be completely 
reduced to the one considered here since it is inherently two 
dimensional (see Section f3. 31 and Appendix [Cj. The results 
can therefore be expected to be quantitatively different, as 
one can readily verify. Indeed, along a magnetic surface of 
a steady jet the energy flux per unit rest mass flux is con- 
served, which for a cold flow leads to the integral 

w = r(i + o) , 

where pj = Tj(l + Oj) and where Tj and <jj represent the 
initial jet Lorentz factor and magnetization, respectively. 
Thus, the Lorentz factor at the tail of the rarefaction wave, 
where the magnetization has decreased to zero, has to be 
T = r 3 (l + <jj). In contrast, equations Q and (I10|l imply 
that F = Fj(l + 2a j), where Oj is the magnetization of 
the undisturbed left state of the Riemann problem (which 
we identify with the initial magnetization in the steady-jet 
problem). 



3.1.2 Case 2 

Here we consider a more complicated Riemann problem 
where the velocity of the left state is not zero but has a 
component tangent to the discontinuity, v = (0,0,Vj). This 
corresponds to the jet expansion as seen in the jet source 
frame. The magnetic field of the left state is B — (0, Fjbj, 0), 
where bj corresponds to the jet azimuthal magnetic field as 
measured in the comoving frame. In addition, we no longer 
assume that the left state is cold, which will allow us to 
compare the purely hydrodynamic and MHD cases. In the 
rarefaction wave v = (v x ,0,Vz) and B = (0 , Ft>, 0). This 
case has been analysed in iRomero et. "ail (|2005h . Integration 
of the simple wave equations ([5]) leads to the integrals 



Is = s , 
lb = b/p, 

I z = hgU z 

and 

2 U - v x /c 



(ii) 

(12) 
(13) 

(14) 



where s is the specific entropy, u z = F(v z /c) is the z 
component of the 4- velocity, 

W + 6 2 /47T 



h g = 



pc A 



is the generalized specific enthalpy and 



X(P) 



1+ui 



1+ ut 1 



1/2 



££ 
C 



(15) 



(16) 



Here Cf is the fast-magnetosonic speed in the direction nor- 
mal to the magnetic field as measured in the fluid frame, 
with 



(17) 
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c a is the Alfven speed, with 



c / b 2 + 4nw ' 
c s is the sound speed, with 

1 f dp 



h \dp 



(18) 



(19) 



and h = w/pc 2 is the specific enhalpy. The function X(p) 
is a function of density only and the other variables are 
eliminated via the integrals (|ll|l . (|12[) and (|13|l for a given 
equation of state. The values of these integrals are dictated 
by the left state of the problem. In particular, 

/+ = and I z — fj,j(vj/c), 

where 

Pi = hjTj(l + (Tj) 
and 

a j = bj/4nwj. 

In these expressions we use the index l j' to indicate the left 
state variables as this state corresponds to the undisturbed 
jet. Also, (jj is the magnetization parameter and pj is the 
energy integral that inc orporates a thermal contribution (see 
iKomissarov et aT1l2009l ). 

In this case the characteristics equation (0 yields 



£ - v x 



= -ct [l + t£(l -cf/c 2 )] 



-1/2 



(20) 



This result has a straightforward interpretation. On the left- 
hand side one immediately recognizes the relativistic expres- 
sion for the addition of velocities, whereas the right-hand 
side gives the speed of a fast-magnetosonic wave as mea- 
sured in a frame that moves along the x axis with the same 
speed as the fluid, i.e. v x (see Appendix lA")l , 

The procedure for constructing the self-similar solutions 
using the above analytical results is the following. ^From 
equations (|11[) - (|14[1 and the equation of state, one finds the 
functions w(p), b(p), u z (p), v x (p) and the thermal pressure 
p(p) Finally, equation (|2Up allows one to obtain x/t = £(p). 
In particular we find that the head of the rarefaction wave 
propagets with the speed 

Choad = -C tj [l + U z . (l - c| /C 2 )] , 

whereas the wave's tail advances with the speed 

Ctail = V x (0) . 

Although in general the solutions can only be found 
numerically, one can derive fully analytic results in the ultra- 
relativistic limit, as follows. From, equation 131) one has 



u z = p,j — — . 

C kg 

For highly relativistic jets v z 
more, p/pc 2 — > and b 2 /pc 2 
then yields 



(21) 

i. Vj ~ c and v x -C c. Further- 
as p/pj — > 0. Equation (f2T|) 



4 For a polytropic equation of state with index 7, the enthalpy 
is w = pc 2 + [7/(7 — l)]p, in which case equation becomes 
p/p 1 = const. 



P/Pi 



0. 



(22) 



The Lorentz factor at the tail of the rarefaction wave (i.e. 
at the boundary with vacuum) is thus found to equal the 
value expected in a steady jet, as discussed in Section [3TTTT] 
The generalized one-dimensional Riemann problem consid- 
ered here (Case 2) therefore provides a better representation 
of the simulated two-dimensional jet problem than the sim- 
pler model considered above (Case 1). 

In our jet simulations we deal with cold (c s = 0), super- 
fast-magnetosonic (r^ 3> ^j^ 2 ), relativistic (vj ~ c) jets. In 
this limit we have 



ho = 1 + a. 



h = l, i = 



Equation (I21|l then gives 
r - 

1 + a' 

whereas equation ()12p yields 
P 

a = a-j — . 

Pi 

Combining the last two results we find that 



(23) 

1 + CTjp/pj 

Moreover, in this limit the integral in equation (|14[) assumes 
a simple analytic form. Specifically, when c s = and v z ~ c 
we have u z fa T = p,j/(l + cr) and thus 

' r<l_L \3, 211/2 



X{p) 



(l + a) 2 + p 2 



In the super-fast-magnetosonic regime Tj 3> u^' 2 and p 2 — 
(1 + o-j) 2 F 2 > (1 + aj) 3 . This makes it possible to simplify 
the expression for X(p) even further, 



Hp) 



which leads to 



Pi 



1/2 



In 



\C-Va 



2a 



1/2 



Pi 



1/2' 



This shows that v x -C c, which enables us to approximate 
the left side of this equation as v x /c and write 



c 



1 2a 



1/2 



r, 1 + a* 



i-l^ 

Po 



1/2 



(24) 



For the jet problem this implies that the jet opening half- 
angle 8 V increases only by the amount 



A6 V 



1 2a 



1/2 



r,- 1 + a, 



(25) 



These considerations imply that the product F# v can in- 
crease to values that are significantly greater than 1 and 
that this occurs mainly on account of the increase in V. 
This is indeed the behaviour seen in our jet simulations. 

3.2 Hydrodynamic versus magnetic mechanisms 

The rarefaction acceleration process has been seen in numer- 
ical simulations of both unmagnetized and magnetized flows 
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Figure 4. Simple rarefaction wave solutions of rclativistic MHD. Each column corresponds to a particular Riemann problem. The 
parameters of the left state are given at the top of the column and the right state is always vacuum. In all cases fij = 600 and the polytropic 
index is 7 = 4/3. The last row shows the Minkowski diagrams for the solutions. Here the colour image represents the distribution of the 
Lorcntz factor (color) and the contours show the worldlines of fluid parcels initially located at x-i = —1, —0.8, —0.6, —0.4, —0.2, —0.02, 0. 
The head of the rarefaction wave can be seen on these plots as the location where the Lorentz factor starts to grow and the worldlines 
start to bend to the right. 



and a number of conclu sions have been draw n about the role 
of the magnetic field (jMizuno et al.ll200ct ). As both cases 
may have applications in astrophysics, we are motivated to 
extend our analysis and investigate the role of the magnetic 
field in the rarefaction mechanism a bit further. In particu- 



lar, it is interesting to see if the presence of strong magnetic 
fields can lead to some observationally identifiable features. 
To check on ths we have derived self-similar solutions of the 
Riemann problem described in Section [3. 1.21 for four differ- 
ent left states: a cold MHD flow, a hot MHD flow and two 
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purely hydro-dynamic (HD) flows with different values of Tj. 
In all of these cases the magnitude of the energy flux per unit 
rest-mass energy flux pj is kept the same. The Lorentz factor 
near the boundary with the vacuum approaches this value, 
independently of the other characteristics of the rarefaction 
wave (see equation I22p . However, the spatial distribution of 
the Lorentz factor near the right boundary turns out to be 
quite sensitive to the magnetization of the left state. 

Fig. [4] shows the results of the numerical integrations 
for an ideal gas with poly tropic index 7 = 4/3. If we com- 
pare the properties of the cold MHD flow (first column) with 
those of the HD flow (third column) we see that the accel- 
eration is much stronger in the cold MHD case. In the first 
place this is due to the fact that, in the MHD case, the accel- 
eration occurs over a larger volume of the jet. This is simply 
a reflection of the difference in the speeds of the rarefaction 
wave's head - in the HD case it is determined by the sound 
speed c s < c/v3, whereas in the MHD case it is determined 
by the fast-magnetosonic speed Cf, which can approach the 
speed of light. 

A second reason for the difference between the MHD 
and HD cases has to do with the dependence of the (gen- 
eralized) specific enthalpy on the mass density. Whereas 
b 2 /Anpc 2 is proportional to the mass density p, w/pc 2 is pro- 
portional to p 7-1 = p 1 ^ 3 . As the density drops across the 
rarefaction wave, the magnetic energy declines much faster 
than the thermal energy, which explains the relative ineffi- 
ciency of the hydrodynamic case. For example, in the HD 
flow shown in the third column of Fig. SI the density is 
p = pj/125 at x/ct = 0.01, but the value of w/pc 2 is still 2 
(meaning that the Lorentz factor is only a half of its max- 
imum value, F = Pj/2). In contrast, in the cold MHD case 
b 2 /iivpc 2 is only 1/25 at the same value of the density, mean- 
ing that r is already ~ pj. Note that the time dependence 
of the density of a particular fluid parcel is different in the 
two cases. As shown in Appendix [B] in the MHD case the 
density drops with time as p oc T~ 2 ^ 3 , whereas in the HD 
case it declines faster (p oc J 1-6 / 7 ) due to the higher speed 
of the tail of the rarefaction. Nevertheless, the (generalized) 
specific enthalpy in the MHD and HD cases is proportional 

0/3 1/3 — 2/7 

to p oc T ' and p 1 ccT , respectively, implying that 
the increase of the Lorentz factor is much stronger in the 
former case. 



3.3 The role of the bunching function 

The above interpretation of the rarefaction-acceleration 
mechanism is consistent with the general analysis of 
the acceleration of co ld, steady-state jets presented in 
iKomissarov et a l. (2009), where it was shown that the bulk 
flow acceleration is intimately related to the shape of the 
field (or stream) lines through the 'bunching function' 



p-F 

r 3 



d(Tv) _ 2 dS 
dt mc di ' 



(27) 



7vr 2 B p 



(26) 



where r is the cylindrical distance from the rotation axis, 
Bp is the amplitude of the poloidal magnetic field and 'J = 
J Bp ■ dS is the magnetic flux function. In the super-fast- 
magnetosonic regime the Lorentz factor increases when S 
decreases. Specifically, one can show that 



where I is the distance measured along the poloidal mag- 
netic field lines and the effective rest-mass m is constant 
on magnetic flux surfaces (see Appendix [C}. The coeffi- 
cient in front of the derivative on the left hand side of this 
equation vanishes at the fast-magnetosonic critical surface 
which, in the highly relativistic limit, implies F ~ p 1 ^ 3 - 
Thus, in the super-fast-magnetosonic regime, acceleration 
corresponds to a decrease in S, whereas in the sub-fast- 
magnetosonic regime it corresponds to an increase in the 
bunching function. This is analogous to the transonic hy- 
drodynamic flow in a de Laval nozzle, with 1 /S playing the 
role of the nozzle cross section. 

Equation ([26} shows that, for S to decrease, B p should 
decrease faster than r -2 , i.e. the separation between neigh- 
bouring magnetic flux surfaces should increase with distance 
faster than r, their cylindrical radius. In c onfined flows, such 
as mo del Bl and the other cases studied in lKomissarov et alj 
(|2009h . this is realised through the stronger collimation of 
the inner flux surfaces relative to the outer ones. For this 
reason, the acceleration mechanism at work in such flows 
can be dubbed the collimation mechanism. 

The acceleration mechanism that operates in models 
Bla and Bib during the transition to the ballistic regime is 
different from the collimation mechanism in that it involves 
a rarefaction wave, and we therefore dub it the rarefac- 
tion mechanism. In this process, the rarefaction wave that is 
launched at the jet boundary at the point where the chan- 
nel widens reaches the jet axis much further downstream. 
Therefore, the outer field lines straighten much closer to the 
source than the inner ones (see Fig. [SJ. The net effect is 
again that Bp decreases faster than r~ 2 . Inspection of mod- 
els Bla and Bib shows that the magnetic bunching function 
indeed decreases along the magnetic field lines (see Fig. 
and that in both of these models the flow becomes super- 
fast-magnetosonic well upstream of the channel widening 
point0 



4 DISCUSSION AND CONCLUSIONS 

Our numerical simulation s confirm the discovery by 
iTchekhovskov et all (|2009bh of the additional fast accelera- 
tion of relativistic jets during the transition from a confined 
to an unconfined regime. We have analyzed the acceleration 
mechanism and concluded that it is related to the sideways 
expansion of the jet that is triggered by the elimination of 
support from the confining wall. This expansion does not 
lead to a significant increase of the jet opening angle be- 
cause of the high Mach number that the flow has already 
attained by the time it reaches the location where the chan- 
nel opens up. The induced rarefaction wave is nevertheless 
strong enough to lead to the conversion of a large fraction of 



J The sub— fast-magneto sonic regime is applicable in model E of 
IKomissarov et al] (120091 ') . for which it was found that, in contrast 
with the super-fast-magnetosonic regime considered in this paper, 
a widening of the jet boundary does not lead to a significant 
acceleration of the flow. 
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Figure 5. Straightening of the jet magnetic field lines after the 
channel opens up. In this sketch the solid lines show the magnetic 
field of an externally confined jet, the dashed lines show the radial 
magnetic field lines of a ballistic jet, and the dotted line shows 
the location of the rarefaction wave front. The magnetic field lines 
of a jet that becomes free after encountering a sudden widening 
of the confining wall (point A) are represented by the solid lines 
upstream of the rarefaction front and by the dashed lines down- 
stream of it (where we assume for simplicity that the field lines 
straighten out right after crossing the front). The figure shows two 
such field lines, which cross the rarefaction front at points A and 
B, respectively. As one moves along the jet from A to B the inner 
line is still parabolic whereas the outer one is already straight. 
Thus, the separation between these lines increases faster than r 
and the amplitude of poloidal magnetic field decreases faster than 
r" 2 . 



the remaining magnetic energy into kinetic energy of bulk 
motion. 

The rarefaction-acceleration mechanism and its po- 
tential relevance to relativ istic astrophysica l flows were 
previously considere d by lAlov fc Rezzollal (|2006l ) and 
M izuiio et all |200Sl ) for unmagnetized and magnetized 
flows, respectively. As discussed in Section [331 the slow ex- 
pansion of the left front of the rarefaction and the weak de- 
pendence of the enthalpy on the density in a hydrodynamic 
outflow, as compared to the faster expansion and stronger 
dependence of the magnetic energy on the density in MHD 
cases, is evidently behind the finding of lMizuno et all (|2008l ) 
that this acceleration mechanism is more efficient in MHD 
outflows than in purely hydrodynamic ones. 

Although this mechanism operates even in Newtonian 
flows, its role there is not that important. This is because 
a high fast-magnetosonic Mach number in the Newtonian 
regime implies that most of the magnetic energy has already 
been converted into kinetic energy of bulk motion, so the 
effect of any additional acceleration is rather insignificant 
even if it involves a full conversion of the remaining magnetic 
energy. The situation is different in the relativistic case since 




W I , i , i , i 

° 2 4 6 

lo Sio z /r lc 

Figure 6. Evolution of the bunching function 5 along the mag- 
netic surface with ^ = 15 for models Bl (solid line), Bla (dash- 
dotted line) and Bib (dashed line). 

the flow can remain magnetically dominated even in the 
highly super-fast-magnetosonic regime. 

Note in this connection that the gradual increase of the 
Lorentz factor of an unmagnetized long/soft GRB outflow 
as it first emerges from u nderneath the surface of the associ- 
ated progenitor star (e.g. iMorsonv et al.ll2007r ) has a differ- 
ent origin than the rarefaction acceleration considered here. 
In the unmagnetized case, a 'cocoon' of moderately rela- 
tivistic, shocked jet and stellar material that formed while 
the outflow was confined within the star is the first compo- 
nent to be revealed. The second component to emerge is the 
faster shocked jet material that occupies the region between 
the jet head and the reverse shock that was induced by the 
interaction of the outflow with the stellar envelope. Finally, 
the unshocked, high- velocity jet itself comes out. The be- 
haviour of the latter component after it leaves the star is 
initially influenced by the presence of the 'shielding' cocoon 
gas. Whether significant rarefaction acceleration could oc- 
cur under these circumstances even in a jet that reaches 
the stellar surface with a measurable fraction of its thermal 
energy not yet converted into kinetic energy remains to be 
investigated. 

We have shown that when the rarefaction wave propa- 
gates inside the jet it produces a region where the separation 
between neighbouring poloidal magnetic surfaces increases 
faster than their cylindrical radius. This is a common fea- 
ture of magnetic acceleration in the super-fast-magnetosonic 
regime and is also a characteristic of the magnetic accelera- 
tion process of externally confined flows by the collimation 
mechanism. However, there are also significant differences 
between these two mechanisms. Most importantly, acceler- 
ation by the collimation mechanism can be sustained over 
many decades in distance and produces asymptotic flows 
with a ~ 1, whereas the rarefaction mechanism continues 
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only for as long as it takes for the rarefaction wave to cross 
the jet (of the order of the current jet length in our numer- 
ical models), and subsequently the flow enters a phase of 
ballistic expansion wherein the magnetic acceleration is so 
slow that it can be ignored. 

We emphasize that the present study has involved 
axisymmetric simulations. A fully three-dimensional (3D) 
treatment is needed to address the issue of stabil- 
ity, although we note tha t the results to date (e.g. 
iMcKinnev fc Blandford1l2009h indicate that magnetically ac- 
celerated relativistic jets basically remain stable as they 
propagate away from the source. It is also not yet entirely 
clear whether the acceleration/collimation properties of 3D 
outflows remain the same as in the axisymmetric case, but, 
at le ast on the basis of Newtonian simulations (e.g. iMollI 
I2009D . it appears that the acceleration efficiency remains 
roughly the same. 

The most effective magnetic acceleration of steady flows 
could be achieved via a combination of both mechanisms - 
first the collimation mechanism produces a flow with a ~ 1, 
and then the rarefaction mechanism provides additional ac- 
celeration, resulting in a particle-dominated {a < 1) flow. 
Model Bib is an illustrative example of such a combina- 
tion (see Fig. [5J, and it is interesting to examine whether 
it could in principle be realized in long/soft GRB sources. 
If the jet originates from a rapidly rotating black hole then 
the light cylinder radius is of the order of a few gravitational 
radii, r\ c ~ 2x 10 6 cm for a 3 M© black hole. If the confining 
medium is the stellar envelope, the collimation mechanism 
operates up to ~ 10 5 n c . A quick inspection of Fig. [5] shows 
that this corresponds to the channel widening point in our 
model Bib. Thus, it may indeed be possible to accelerate 
long/soft GRB flows to the high Lorentz factors indicated by 
observations, with a decreasing to values < 1 (as happens in 
this example on scales z ~ 10 12 cm). However, it seems un- 
likely that the low magnetization (a -C 1) required for effec- 
tive dissipation in the internal-shocks model for the prompt 
7-ray emission can be attained in this scenario, which sug- 
gests that magnetic energy dissipatio n might need to be in- 
voked to explain this emission (e.g. iLvutikov fc Blandfordl 
I2003T ). If the initial magnetization cjq is very high, well above 
~ 10 3 , the jet will still be magnetically dominated when it 
enters the ballistic regime. 

As pointed out by iTchekhovskov et al.1 (|2009bf ). the 
rarefaction-acceleration mechanism can in principle give rise 
to magnetically accelerated outflows that satisfy T6j 2> 1. 
In the standard model of uniform GRB jets, this condi- 
tion allows a late' (on a time scale > 1 day) panchro- 
matic jet break in the afterglow light curve. However, as 
the asymptotic structure of magnetically accelerated jets is 
far from being uniform, further investigation is required to 
fully establish this result. Furthermore, in some of the GRB 
sources detected by Swift there have been indic ations of an 
early (on a time sca l e of > 1 hr) jet break (e.g. iPanaitescul 
120071 : iKamble et ail I2009T K corresponding to T6j < 1 and 
9j < 1°, as predicted by pure collimation-acceleration mod- 
els {k omissarov et al.l l2009h . If these inferences are cor- 
roborated by further observations, the reasons why certain 
long/soft GRB sources show no evidence of rarefaction ac- 
celeration will need to be clarified. It is, however, also con- 
ceivable that the acceleration of jets that exhibit a late break 
in the afterglow lightcurve is predominantly thermal, since 



in this case the value of Fdj is not limited as in the mag- 
netic acceleration scenario. One may, however, be able to 
distinguish between the thermal and magnetic mechanisms 
- at least in some cases - based on other observational diag- 
nostics, such as the appearance o f a photospheric emission 
component (e.g. iRvde et al.ll2010l ). This issue could poten- 
tially be further illuminated by studies of short/hard GRB 
afterglows. Since the outflows in the latter sources evidently 
do not propagate through a stellar envelope, the rarefaction- 
acceleration mechanism would not operate in this case (or 
at least not in the same way as in long/soft GRB sources). 
If short/hard GRB outflows are accelerated predominantly 
by magnetic stresses, this suggests that there should be no 
late-time breaks in their afterglow lightcurves (or at least 
fewer than in the case of long/soft GRBs). While th e data 
on sh ort /hard GRB afterglows are still sparse (e.g. iNakarl 
I2007I) . this prediction could in principle be tested as more 
such afterglows are observed. 
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APPENDIX A: INTERPRETATION OF 
EQUATION pOj) 

The derivation of the characteristics equations that yield 
the phase speeds of relativisti c MHD waves can be found 
in a n umber of referenc es (e.g. lAnile fc Pen nisi 1987; lAnilel 
1 19891 ; iKomissarovl I1999T ). For the fast-magnetosonic waves 
the result is 



1-vi 



r 4 (A-^) 4 (l-c 2 ) + (l-A 2 ) x 
[{<2/w g )(V - A6 ) 2 - T 2 (A - v x ) 2 c 2 ] = , 



(Al) 



where w g = w + b 2 /Atv is the generalized enthalpy and b" 
is the magnetic fi eld 4-vector (see equations 29 and A20 in 
IKomissarovl 1 19991 ) . In order to have compact equations we 
employ in this appendix geometric units in which c = 1. The 
components of b" can be obtained given the magnetic field 
and flow velocity 3-vectors via 



b° = TB i v i , b 1 = (B 1 -Vv l b°)/Y. 



(A2) 



In the problem under consideration b° = and b x — 0. 
Moreover, in the frame comoving with the fluid along the x 
axis v x = and equation (|A1|) yields 



A 2 = c 2 [f 2 + c 2 (l — f 2 )] 1 



(A3) 



Here we use a tilde to indicate quantities in the comoving 
frame (comoving only in the x direction). Using the equa- 
tions of Lorentz transformation we find 



r = rr(i - vv x 



(A4) 



where Vi and F are the flow parameters as measured in the 
frame of the initial discontinuity and where F and v describe 
the motion of the comoving frame. Since v x = v, this equa- 
tion yields 



1-vl-vi 



1 + ui, 



(A5) 



where u z = Yv z . Substituting this result into equation (|A3|) . 
we obtain 



A 2 = c 2 [l + u 2 (l - c 2 )] 1 



(A6) 



Thus, the right-hand side of equation Q20p is indeed the fast- 
magnetosonic wave speed as measured in a frame that co- 
moves with the fluid along the x axis. 



APPENDIX B: THE RAREFACTION WAVE IN 
THE ULTRA- RELATIVISTIC LIMIT 



The Riemann problem considered in Section 13.1.21 is de- 
scribed by equations (|11|| - (|14|) and (|20|l . Here we add that 
one can follow the wordline of a fluid parcel initially located 
at Xi(< 0), by integrating the equation 



dint 



(Bl) 



which follows from the definition £ = x/t. The dimensionless 
time T = ct/(—Xi), for times T > T; = — l/£head (i.e. after 
the head of the rarefaction wave has already passed Xi), is 
given by 



T = Ti exp 



Pi 



d p 



dp . 



(B2) 



Bl The ultra-relativistic cold MHD limit 

The expressions for the Lorentz factor V and the expansion 
speed v x as functions of the density p in the cold (c s = 0), 
ultra-relativistic (v z ~ c) limit have already been derived 
in Section 13.1.21 and are given by equations (|23p and Q24p . 
respectively. Using the same simplifications, we can find the 
density and the Lorentz factor as functions of £ = x/t. Equa- 
tion H20|) can be writen as 



x 
ct 



1 

l< j 



2a 



1/2 



Pj 



1/2 



Pj 



3/2 



(B3) 



' /Tj for the head and 



implying (x/ct) min = 

(x/cf) ma x = 2o "j /2 

The Lorentz factor is V — /ij / (1 + (Tjp/pj), with 



2a- /jij for the tail of the rarefaction wave. 



■ sinh 



- arcsinh 



(of 



El — 

2 ct 



(B4) 



found by inverting equation (|B3[) . We can also express the 
Lorentz factor as a function of the dimensionless time T, 
since equation ()B2|| implies 



(B5) 



-l/3i-i2/3rrr-2/3 

The resulting expression is 
Mi 



r = 



l+o- 



2/3 r 2/3 



p^/oj 1 _ 2 /3 



(B6) 
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B2 The ultra-relativistic HD limit 

Similar approximations in the purely hydrodynamic limit 
{a j = 0) imply 



fi.i 



1 + Qn - 1) (P 



c 



1 [(2 - 7 ) (hi — 1)] — I [(2 - 7 ) (h 3 - 1) g-r- 1 ] 
W (7-l) 1/2 (2-7) 1/2 

I [C] = C 1/2 (i + C) 1/2 + in [C 1/2 + (l + C) 1/2 ] 



c7 



c 



,7-1 



-(7+l)/2 



( 7 - 1) (hj - 1) Q 

l + (2- 7 )(ft J -l)^- 1 

l + (2 - 7 ) (ftj - l) e 1 - 1 



(7-1)^-1) 



1/2 



1/2 



(B7) 

(B8) 
(B9) 

,(B10) 



where q = p/pj. (Note that, for 7 — > 2 and hj — > 1 + <jj, we 
recover the cold MHD case considered in Section [B2]) 

For 7 = 4/3 and T significantly larger than OApj, equa- 
tion (|B10[) simplifies to 



^-3/rp6/7y-6/7 



(Bll) 



Subtituting this result into equation ()B7[) . we infer the de- 
pendence of the Lorentz factor on T, 



r = 



/'./ 



i + 3 1 / 7 (/i J -i) 6/7 r 2/7 T- 2 / 7 



(B12) 



Comparing the results (|B6|) and (|B12|l . we see that the 
rarefaction-acceleration process is faster in the MHD case 
than in a purely hydrdynamic flow. This is part of the reason 
for why this mechanism is more efficient in MHD jets, as 
discussed in Section T3. 2 1 



APPENDIX C: THE BUNCHING FUNCTION 
AND MAGNETIC ACCELERATION 

At distances from the central source where the flow can be 
considered cold and the azimuthal velocity small, the com- 
ponent of the momentum equation along the motion can be 
written as 



^(i + r^ 



di 



B$ d(rB^ 
4irr di 



(CI) 



where p is the rest-mass density, r is the cylindrical coordi- 
nate, is the azimuthal component of the magnetic field, 
I is the arclength along a poloidal field (or stream) line, and 
derivatives with respect to £ are taken along a give n field 
line (see e.g. equation 20 in lVlahakis fc Konig]||200 3 ). 

We can further simplify equation (|C1|) by using two in- 
tegrals of motion, the mass flux per unit magnetic flux (k) 
and the angul ar velocity of magnetic fi eld lines ($1) (see e.g. 
Section 2.1 in I Komissarov et alj|2009h . Using their expres- 
sions in the limit of a small azimuthal velocity, we have 



k = 



pYv 

Br, 



r B v 



(C2) 



Equation (|C1[) can be written in the form of a momentum 
equation of a point particle in a potential field: 



d(Tv) 
1 dt 



dV 

"dl 



(C3) 



The effective rest mass, m, is inversely proportional to the 
magnetization, 



S-K 2 kc 
*ft 2 



(C4) 



and the effective potential energy is related to the bunching 
function, 



V = 



25 

v/c 



(C5) 



The corresponding energy equation can be written as 

mrc 2 + V = rape 2 , (C6) 

where the integral of motion p, represents the total energy 
flux per unit rest-mass energy flux. 

The bunching function S is directly connected to the 
geometry of the flow. The cross sectional area between two 
neighbouring flux surfaces and 9 +6^ is <5\l/ / B p . Thus, the 
function S is proportional to r 2 over this area and decreases 
whenever the flow expa nds in such a way that this area 
increases faster than r 2 dVlahakidl2004 see also Section 5.1 
in lKomissarov et al.ll2009t ). The evolution of S is determined 
by the transfield component of the momentum equation; it 
depends on the external pressure that confines the jet (or, 
equivalently, on the shape of the 'wall' that defines its outer 
boundary). 

Since the potential depends not only on S but also on 
the velocity, it is useful to separate out the effect of geometry 
in the equation of motion. Combining equations (|C3|l . (|C5|l 
and (IC6I) wc can write 



H-T\ d(Fv) 



dC 



2 dS 

mc dl 



(C7) 



The resulting critical point corresponds to the fast magne- 
tosonic surface, where dS/di = and V ~ p}^' A (in the 
relativistic regime where 1 -C Y -C T 3 ). Acceleration in 
the super-fast (sub-fast) magnetosonic regime corresponds 
to decreasing (increasing) 5, respectively. 

If at some point the curvature of the jet boundary sud- 
dently increases, the adjustment of the magnetic field from 
the old to the new curvature corresponds to a fast expansion 
and decrease of the bunching function. This is precisely the 
effect of the rarefaction wave analysed in Section[3] As an ex- 
ample, suppose that the initial shape is parabolic, 2 = c\r C2 , 
and the final shape is conical, 2 = zq + r/tani9 tr , as shown 
in Fig. [S] The variables 20 and #tr are functions of ^ and 
can be found from a smooth matching of the magnetic field 
along the transition surface 2 tr = 2t r (r , tr ) (dotted line in 
Fig. [SJ. Downstream from this surface the bunching func- 
tion 5 declines as 



S tI + AS 



1 + 



AS rt r 
<Str r 



AS ■ 



-Su- 



shi #tr dZo 
rtr d$tv 



(C8) 



(see equation 35 and the related discussion in lVlahakisll2004l 
for further details). The appearance of the variable 20 (*T/) 
in the equation describing the conical boundary (2 = 20 + 
rj tani9 t r) is crucial for enablig the additional acceleration 
(by making AS < 0): this is elaborated in the discussion 
of the difference between 'type I conical' (in whic h zp = 0) 
and 't ype la conical' (in which 20 7^ 0) shapes in IVlahakia 
(|2004h ; similar categories exist for parabolic shapes as well 
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(type I, la, II and Ha). Although the adopted channel shape 
in model Bib is a simple cone whose vertex is located at the 
origin, the field lines clearly have a 'type I conical' shape, 
i.e. their projection crosses the r = axis at zo < with 
dzo/d^ > (see the right panel of Fig. [1] and Fig. [5). 

Equation (|C8|I fits well the curves seen in Fig. [5] In 
one decade or so in cylindrical distance downstream from 
the transition radius r tr , the bunching function drops from 
Str to Sn + AS. The difference AS (< 0) corresponds to a 
difference AV ~ 2A<S in the potential energy and thus to an 
acceleration Ar = — 2AS /mc 2 . 
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